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ABSTRACT 



Polycyclic aromatic hydrocarbons (PAHs) produce characteristic infrared emission bands that have been observed in a wide range 
of astrophysical environments, where carbonaceous material is subjected to ultraviolet (UV) radiation. Although PAHs are expected 
to form in carbon-rich AGB stars, they have up to now only been observed in binary systems where a hot companion provides a 
hard radiation field. In this letter, we present low-resolution infrared spectra of four S-type AGB stars, selected from a sample of 
90 S-type AGB stars observed with the infrared spectrograph aboard the Spitzer satellite. The spectra of these four stars show the 
typical infrared features of PAH molecules. We confirm the correlation between the temperature of the central star and the centroid 
wavelength of the 7.9 /im feature, present in a wide variety of stars spanning a temperature range from 3 000 to 12 000 K. Three of 
four sources presented in this paper extend this relation towards lower temperatures. We argue that the mixture of hydrocarbons we 
see in these S-stars has a rich aliphatic component. The fourth star, BZ CMa, deviates from this correlation. Based on the similarity 
with the evolved binary TU Tau, we predict that BZ CMa has a hot companion as well. 

Key words, techniques: spectroscopic - stars: AGB and post-AGB ~ stars: circumstellar matter - stars: S-type stars 



1. Introduction 

Polycyclic aromatic hydrocarbon (PAH) molecules are large 
molecules containing carbon and hydrogen atoms. They carry 
typical infrared (IR) emission featur es that have b een observed 
in many astrophysical environments (Tielensll2008h . Since these 
features are generally attributed to the IR fluorescence of 
ultraviolet-pumped molecules, we expect PAH features to arise 
I from regions where c arbon -rich material is exposed to ultra- 
■ violet (UV) radiation dLege r & Puget 19 841: ICohen etalj|1985l; 
! |Puget& Leger"l989': Allamandola et al.ll989l) . PAHs plav a ma- 
' jor role in photoelectric heating processes and the ionizat ion bal- 
[ ance of the interstellar material dLepp & Dalgarnolll988h . 

Although interstellar PAH molecules are thought to originate 
in the winds of carbon-rich asymptotic giant branch (AGB) stars, 
there is little observational support for this idea. In a complete 
ISO/S WS survey of 50 carbon-rich AGB stars, iBoersma et alJ 
(|2006|) detected PAH emission in only one source, TU Tau, a 
carbon-rich AGB star with a hot A2 companion providing UV 
radiation. Recently, Sloan et al. (2007) has detected PAH emis- 
sion in an R-type carbon-rich giant with a circumstellar disk, 
which is the likely region where the PAH emission originates. 
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This lack of PAH detections in carbon-rich AGB stars is most 
likely the result of the high optical opacities of carbon-rich dust: 
when formed, the PAH molecules are shielded from optical and 
ultraviolet photons by the opaque carbonaceous circumstellar 
material (Jageret al. 1998). 

S-type stars are objects on the ascent of the AGB. They pass 
through the phase where the photosphere turns from oxygen- 
rich to carbon-rich (C/O ~ 1). Because S-type stars are not 
yet forming carbon-rich dust, they have less opaque circumstel - 
lar shells than their carbon-rich successors dJager et al.ll2003h . 
Nevertheless, due to shock-induced non-equilibrium chem- 
istry effects, S-type stars could display carbon-rich molecules 
(Cherchneff 2006). Recent studies of ISO/SWS spectra of S-type 
AGB stars indeed show the presence of m olecules hke HCN and 
C2H2 dYang et al.ll20"07tlHonv et al.ll2009h . 

We have studie d the data obtained with the Spitzer Infrared 
Spectrograph (IRS. iHouck et al.l l2004') of a sample of 90 intrin- 
sicQS-type stars (Program ID 3 073 7, PL S. Hony). A description 
of the sample can be found in lCami et al.l (I2009h . They present 
the detection of SiS absorption bands, in a large subset of the 
sample. In this paper, we present the detection of PAH emission 
around four cool S-type AGB stars. Three of the four objects ex- 
tend the known correlation between the centroid wavelength of 



' Intrinsic S-type stars are not enriched in C or s-process elements by 
accretion from an evolved companion, but through internal nucleosyn- 
thesis 



2 



K. Smolders et al.: When an old star smolders: 



Table 1. Anchor points of the continuum spHnes. 



Feature 


Anchor Points 






Class B/C sources 


Class A sources 


/jm 




fim 


6.2 


5.7,5.9, 6.1,6.5,6.7, 6.9 


idem 


7.9 


7.0, 7.2, 7.4, 9, 9.2, 9.4 6.6, 


6.8, 7.0, 9.0, 9.2, 9.4 


11.2 


10.3, 10.5, 10.7, 11.8, 12.0, 12.2 


idem 



C-C stretch : C-C stretch : C-H out pf plahe bendihg 

: : C-H "in' plane ben'dihg :' ' ' " "C-'H dut of plane bending 



the 7.9 ijm PAH feature and the stellar effective temperature to 
lower temperatures and more redshifted features. We argue that 
the hydrocarbon molecules have a high aliphatic/aromatic ratio. 
Since AGB stars are important producers of the dust in the inter- 
stellar medium, we might be looking at the hydrocarbon mixture, 
as formed in the current day wind. 

2. Observations and data reduction 

The Spitzer IRS spectra were obtained in the low -resolution star- 
ing mode, which places the star in two nod positions. This is first 
done in the short-low (SL) apertures, and then in two long-low 
(LL) apertures. The resulting spectrum covers the 5 to 37 fj.m re- 
gion with a typical resolution of R~100. 

For the data reduction we used the FEPS pipeline, devel- 
oped for the Spitzer legacy program "Formation and evolu- 
ti on of planeta r y syst ems". A detailed description can be found 
in iHines et alj ( l2005 l). For the extraction we used the interme- 
diate droop data products as delivered by the SSC, together 
with the SMART re duction package tools described in detail in 
iHigdon et al.l (|2004|) . The spectra are background-corrected by 
subtracting the images of the two slit positions for each module 
and order The IRS s pectrum of all stars is d e-reddened using the 
extinc tion maps from Schlegel et alj (Il998h and the IR extinction 
law of IChiar & TielensI (120061k 



3. The observed emission features 

The PAH molecules exhibit very characteristic emission features 
in the 5-14 fim region, the most prominent at 6.2, 7.9, 8.6, and 
11.2 fj.m. These features arise from the bending and stretching of 
the carbon and hydrogen bonds in the large molecules. 



3.1. Extraction of the PAH features 

The S-type stars in the sample show a wide variety of infrared 
spectra, containing photospheric absorption bands of e.g. CO, 
H2O, SiO, and SiS. The absorption bands in combination with 
the weak emission features make it difficult to estimate an accu- 
rate continuum. The method of extracting possible emission fea- 
tures combines the use of cubic splines for the 6.2 and 1 1.2 fim 
features and comparison stars for the 7.9-8.6 fim complex. 

For the regions around 6.2 and 1 1.2 yum, the continuum flux 
was estimated by fitting a cubic spline (with smoothness zero) 
through anchor points at different positions next to the features 
(see Table [TJ. By applying random variations of the spectra 
within the calculated error bars and shifting the anchor points 
2 spectral points outwards, we obtained a large set of continua. 
The spread on these continua was used to estimate the error in- 
troduced by fitting a spline to the data. 

Strong water vapor absorption can mimic the prese nce of 
a 6.2 Atm emission feature (see for example Fig. A. 23 in ICamil 
I2OOI . However, the water vapor bands can be easily recognized 
by comparing the full spectrum with synthetic spectra, and are 




Wavelength [^m 



Fig. 1. Top panel: The infrared spectra of the 4 stars (black) and 
the continuum estimates (red) with the characteristic infrared 
emission features from PAH molecules, shown in Rayleigh units. 
As a reference, we show the spectrum of CF Gem, a naked pho- 
tosphere, with estimated continuum. The spectrum of CSS 757 
shows weak PAH emission, as well as a dust emission feature at 
1 1 jjm. Middle panel: The residual flux. Bottom panel: Rescaled 
ISO-SWS spectra of prototypes of the different PAH classes: 
NGC 133 3-3 (Class A), HP 44179 (Class B) and CRL 2 688 
(Class C) (iPeeters et al .1120021: Ivan Diedenhoven et alj|2004l) . 



generally weak in S-type stars. We excluded five stars in which 
the detection of a 6.2 /im feature was due to this effect. 

BZ CMa is the only star in our sample of S-type stars with 
possible PAH emission, which does not show the SiO absorp- 
tion band in the region containing the 7.9 and 8.6 jjm PAH emis- 
sion features. For this star, we used a spline at the anchor points 
mentioned in Table [T]to estimate the continuum. The other stars 
contain the SiO band head, which varies in depth throughout the 
sample of S-type stars. Because this absorption band is asym- 
metrical, a spline is not a good estimate of the photospheric con- 
tinuum. Instead, we used the 40 out of 90 sample stars without 
dust emission (naked photospheres) as template stars, fitted in 
the regions 7.4-7.9 yum and 9-9.4 yum where the contribution of 
the PAH emission is negligible. We used the spread on the 5 
best-fitting template stars as a conservative estimate of the error 
caused by fitting these spectra. 
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Table 2. The 4 S-type stars with significant emission features at 6.2, 7.9, and 11. 2 yum. 





Stellar parameters 


6.2 yum 






7.9 yum 






1 1.2 yum 




Name 




C/O" 


Flux" 


Class"^ 




^cen 


Flux" 




Flux" 






K 




erg cm"^ s"' 






yum 


erg cm-^ s"' 




erg cm"^ s"' 




CSS 173 


3600 ± 150 


0.97+0.01 


3.0 ± 0.7 X 10-' ' 


C 


8.25 + 0.05 


8.38+0.10 


1.6 + 0.2 X 10- 


12 


1.1 +0.2 X 10- 


13 


KRCam 


3400 ± 150 


0.97+0.01 


6.2 + 0.6X 10^" 


C 


8.16 + 0.05 


8.31 +0.06 


1.5+0.1 X 10" 


12 


1.3 + 0.2 X 10- 


13 


BZ CMa 


3400 ± 150 


0.98 ±0.01 


20 + 2 X lO-'-'' 


B 


7.97 + 0.04 


8.02 + 0.04 


3.4 + 0.2 X 10^ 


12 


9.8 + 0.5 X 10- 


13 


CSS 757 


3400 ± 150 


0.97+0.01 


5.6 + 0.7 X 10-13 


C 


8.27 + 0.05 


8.47 + 0.06 


2.2 + 0.2 X 10- 


12 


1.1 +0.3X 10- 


13 



" Derived fro m the comparis on of low-resolution spectroscopy with a dedicated grid of marcs model atmospheres for S-type stars (Plez et al. in 
preparation. iGustafsson et al. 2008) . The low-resolution spectroscopy was obtained for the complete sample of S-type stars, with the William 
Hershel Telescope on La Palma for the sources in the northern sky (range: 4100 A-8300 A; resolution: 2 A) and with the 1.9 m telescope at the 
South African Astronomical Observatory for the southern targets (range: 4900 A-7800 A; resolution: 4 A). The entire dataset will be presented 
in a forthcoming paper. 

" Integrated line flux of the residual flux in the interval [6.12,6.5] fim, [7.4,9.0] yum, or [10.9,11.8] pm. 
According to .Peeters et al.. (,2002) . 



We calculated the integrated flux of the three continuum- 
subtracted features. When all features had a line flux exceeding 
the 3-sigma level, we considered this a positive detection of PAH 
emission. In Table|2]we present the 4 stars with a positive detec- 
tion of PAH emission. Figure[T]shows the spectra, the estimated 
continua and the residual fluxes of these stars. We include the 
spectrum of CF Gem, a star without PAH emission, for compar- 
ison and to indicate the quality of the continuum estimates. 

The detected PAH features in the infrared spectra can have 
an interstellar origin. To exclude this possibility, the background 
spectrum in the off-source nod (i.e. in the other order) was ex- 
tracted. Since no significant features were found in the back- 
ground spectra, we conclude that the PAH emission comes from 
the immediate environment of the star. Moreover, the feature 
profiles are very different from those observed in interstellar 
sources. From here on we focus only on these 4 stars with a 
positive detection of PAHs: BZ CMa, KR Cam, CSS 173, and 
CSS 757. 



3.2. Classification 



iPeeters et al.l (l2002h shows that the C-C stretch modes at 6.2 and 
7.9 jum display a wide spread in peak position. They classify ob- 
jects in three classes from A (bluest features) over B to C (red- 
dest features). Because of the low spectral resolution and signal- 
to-noise of our data, we focused on the broad feature at 7.9 yum. 
For Class C the peak of this feature is shifted to 8.2 fim, mark- 
ing a clear difference with the 7.9 fim feature of Class B PAHs 
and the 7.6 fim feature of Class A PAHs. In Table |2] we show 
the peak wavelengths and the classes we derived for all stars. 
In Fig. [U the residual spectra of the S-type stars with PAH fea- 
tures, i.e. after continuum subtraction, are compared with PAH 
features from each of these 3 classes. 

Based on the peak wavelength, BZ CMa is the only star 
that shows the Class B PAHs, generally found in PNe, Herbig 
Ae/Be stars, and in some red supergiants dSvlvester et al.l[T998t 
[Tielens 2008; Verhoelst et al. 2009). The three other stars are 
members of the small group of known Class C sources, together 
with a few po st-AGB stars dPeeters et al J 120021 : iKraemer et alJ 
l2006tlGielen et al. 2009]), die Herbig Ae/Be source HD 135344B 
(ISlo an et al.' 2005). the intermediate-mass T Tauri star SU 
Aur (Fui'lan et al. 2006), the carbon-rich red giant HD 100764 
dSloanet al.lT2007b . and the oxygen-rich K giant HD 233517 
dJura et alJl2006h . 



4. Correlation between Te^ and feature redshift 

The centroid wavelength (/Icen) is the continuum-subtracted, 
flux-weighted barycenter of the feature. Recently, it has been 
claimed that the centroid wavelength of the 6.2, 7.9, and 
11.2 fim features are shifted to longer wav elengths in targets 
wit h a lower effective stellar temperature dSloan et al.l 120071 : 
iKeller et al.,.2008:.Boersma et al.,.2()08.) . 



4.1. Extending fhe correiation towards cooler stars 

In Fig.|2]we show the centroid wavelength of the 7.9 fim feature 
for all stars in our sample, supplemented by the centroid wave- 
lengths of a comparison sampl e consisting of all Class C sources 
with observed 7.9 iim features ("Kraemer et al. 2006: Furla n et al.l 
2006: Jura et al. 2006: Sloan et al. 2007), the carbon-rich AGB 
star TU Tau (Boersma et al. 2006), the Herbig Ae stars without 
si licate emission from iKeller et al.i (i2008D . and all PAH sources 
in lPeeters etdll d2002h with steUar temperatures below 20 000 K. 
The continuum flux was estimated by fitting a cubic spline 
through the anchor points given in Table [T] as described in Sect. 
13.11 To calculate the centroid wavelengths, we fit a line to the 
left and right of the 7.9 yum feature and removed the 8.6 yum fea- 
ture. Redefining the centroid wavelength for all sources, using 
the same methods as for the S-type stars, allows for a quantative 
description of the correlation. 

Our sample confirms the strong correlation between the ef- 
fective temperature of the central star and the centroid position 
of the 7.9 fim feature. Moreover, the PAH features in the cool 
S-type stars are redshifted the most. The correlation can be ap- 
proximated by a linear relationship for all temperatures below 
12 000 K. A linear regression analysis results in 

/Icen [/^rn] = (8.05+0.01)-(7.54±0.60)x(Teff [10^ K]-0.077) (1) 

with a Pearson correlation coefficient of R - -0.90 and a two- 
tailed probability much lower than 0. 1 percent, indicating a very 
strong correlation. All sources with higher temperatures have 
Peelers Class A features (/Icen = 7.84 + 0.06 yum). A similar cor- 
relation is present for the 6.2 and 11.2//m features, albeit to a 
lesser degree. 

Two stars that do not follow this linear trend are BZ CMa 
and TU Tau. The latter, however, is a known spectroscopic bi- 
nary with a hot A2-type companion (Tefr ~ 9400 K) ( Richer! 
IT971 . If we use this temperature instead of the effective tem- 
perature of the carbon star (2750 K), this star falls within the 
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scatter on the correlation. For BZ CMa we use Eq. ([TJ to predict 
the presence of a hot A7-type companion with a temperature of 
approximately 8 100 K. Additional data is required to check this 
prediction. 

For the three other stars there is no indication of a hot com- 
panion. The detection of PAH features in the IR spectra of these 
stars means that PAHs can be sufficiently excited by stars with 
low effective temperatures, hence soft radiation fields. 

4.2. Aliphatic and aromatic material 



ISloan et aL I (l2007l) propose an explanation for this correlation, 
based on qualitatively similar shifts in laboratory measurements 
of PAHs. Class C PAHs are not subjected to strong UV radi- 
ation fields and thus are relatively unaltered mixes of aliphatic 
and aromatic hydrocarbons. Class A and B PAHs have been 
processed by the hard radiation of the central stars, suppress - 
ing the aliphatic component (e.g. 'Goto et al ."2 0021 12003LI2007I) . 
This agrees with the results of Pino et al. (2008^ who show that 
Class C sources, with a redshifted 6.2 jim feature, have a strong 
a strong aliphatic emission feature at 3.4 yum. 

The three targets following the correlation at very low tem- 
peratures can be classified as extreme Class C sources, with very 
broad, red emission features. CSS 757, KR Cam and CSS 173, 
are thus the objects with the highest aliphatic/aromatic content 
ratio observed in space to date, exhibiting the composition as 
condensed in their winds. 



5. Discussion 

Merely that we can detect hydrocarbon features around cool 
stars without a hot companion is remarkable. These detections 
push the temperature of the radiation field exciting PAHs to 
even lower values than before. Since these stars do not have 
strong UV radiation fields, the hydrocarbon mixture must be 
pumped by optical photons. The theoretical model presented in 
iLi & Draine (2002) predicts that optical pumping is possible for 
ionized PAHs and for large PAH molecules with more than 100 
carbon atoms. Later refinements of this PAH model, based on 
laboratory measurements, show th at even near- infrared photons 
can excite ionized PAH molecules (iMattioda et a l. 2005a b). The 
detection of hydrocarbon emission around the S-type stars im- 
plies that the molecules are either ionized or very large, or, alter- 
natively, that they can be excited by optical photons due to the 
aliphatic component. 

The detection of hydrocarbon emission in combination with 
the SiO absorption band in the stellar photosphere points to 
mixed chemistry. This is unexpected if the ongoing chemistry in 
the circu mstellar environrnent of S-type stars is in thermal equi- 
librium ( iFerrarotti & Gailll2002l) . Because the presence of SiO 
in the photosphere points to an oxygen-dominated chemistry, 
th e majority of C a toms are bound in CO molecules. However, 
in ICherchneffI (l2006i) . the nonequilibrium chemistry in the inner 
wind of AGB s tars is studied, and the re sults differ strongly from 
those found bv IFerrarotti & Gail (12002 ). For stars with an equal 
amount of carbon and oxygen (C/O = 1 ) nonequilibrium, shock- 
driven chemistry predicts molecules in the outer regions of 
the stellar outflo w. These carbon-rich molecules ar e the building 
blocks of PAHs (ICherchneff et al.lll992l;ICaiill2002l) . This is con- 
sistent with all C/O ratios found for these stars being very close 
to unity (see Table |2|l. When compared to our entire sample of 
90 S-type stars, the C/O ratios of these four stars are in the top 
10%, with the others being less carbon-rich. 



This paper 
giants 
post-AG B 
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Fig. 2. The centroid wavelength of the 7.9 ;um feature against the 
stellar effective temperature. BZ CMa and TU Tau are shown as 
black circles. The arrow indicates the shift of TU Tau to the tem- 
perature of the A2 companion. The gray line and area represent 
the correlation, quantified in Eq. ([T]l with the corresponding 99% 
confidence interval. 



6. Conclusions 

In this paper we present a positive detection and identification 
of PAH emission in 4 S-type AGB stars. In this small sample 
we see a clear difference between the strong PAH emission in 
BZ CMa, which can be classified as Peeters Class B, and the 
much weaker PAH features seen in the other sources of Class C. 
We predict that BZ CMa is a binary system with a hot, late-A- 
type companion. 

Our data are consistent with the strong correlation found be- 
tween the centroid wavelength of the 7.9 feature and the 
temperature of the central star. They extend this correlation to- 
wards lower temperatures and more redshifted features. This is 
consistent with the hypothesis that Class C PAHs are hydrocar- 
bon molecules with a high aliphatic/aromatic content ratio found 
around stars with weak UV radiation fields. The hydrocarbons 
around CSS 757, KR Cam, and CSS 173 thus represent the com- 
position as condensed in the AGB wind, before entering the in- 
terstellar medium where harsh UV radiation alters their chemical 
structure. 
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